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At age two, humans are just learning to speak, and they see the 
world in an unbiased way. At age two, the COMPLETE (Coordinated 
Molecular Probe Line Extinction Thermal Emission ) Survey of Star- 
Forming Regions is just beginning to hint at what it will do for our 
world view, and it is unbiased as well. The coordinated observations 
that comprise COMPLETE have already shown us: 1) that large (~ 10 
pc) scale bubbles created by winds from evolved stars dominate maps of 
extended dust emission even in relatively low-mass star-forming regions 
(e.g. Perseus and Ophiuchus); 2) that "giant" (> 1 pc-scale) outflows 
from young stars can be easily detected by applying statistical tools to an 
unbiased CO survey (e.g. we have doubled the number of known outflows 
in Perseus using the Spectral Correlation Function); and 3) that many 
apparent high column-density "clouds" seen in projection are actually 
the superposition of many velocity features (likely at different distances) 
along the line of sight. The upcoming detailed analysis of these new ob- 
servations will quantify the degree to which the apparently ubiquitous 
spherical and collimated outflows can sculpt and stir molecular clouds. 
Taken together with recent numerical simulations' results and measure- 
ments of high velocities for some young stars, the COMPLETE Survey 
also raises, and should soon answer, questions about how, when, and how 
far stars move away from their birthplaces as functions of time. 



1. Introduction 



The prevailing theoretical picture of star formation envisions stars forming inside 
of dense cores, which are in turn embedded in larger, slightly lower-density 
structures. Each forming star is surrounded by a disk, and, when it is very 
young, the star-disk system produces a collimated bipolar flow, in a direction 
perpendicular to the disk. In its broad outlines, this paradigm is very likely to be 
right. In detail, though, many questions concerning the timing of this series of 
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events remain. For example, how long does a star stay with its natal core? How 
long does it remain associated with the lower-density structure (e.g. a filament in 
a dark cloud) where it originally formed? What kind of environment does a star- 
disk system need to keep accreting, or to produce an outflow-and when might 
that reservoir no longer be available to the system? Does a bipolar outflow have 
any effect on star formation nearby? What causes fragmentation into binaries or 
higher-order systems? How much influence do spherical winds (e.g. SNe, B-star 
winds) from previous generations of stars have on the timing of star formation? 
How often is a star in the process of forming likely to encounter an external 
gravitational potential (e.g. from another forming star) strong enough to alter 
its formation process? The complicating issue underlying all of these questions 
is that it is hard to define and understand the long-lasting properties of the 
reservoir from which a star forms if the reservoir itself is highly dynamic. 

The COMPLETE Survey of Star-Forming Regions was born out of a de- 
sire to create an unprecedented comprehensive statistical database with which 
one might have real hope of answering many of these questions. We were in- 
spired to begin COMPLETE at the 2001 Santa Cruz Star Formation Workshop, 
while grousing at a coffee break about how virtually none of the beautiful mea- 
surements of cloud's velocity fields (from molecular line observations), density 
profiles (from extinction measurements), temperature and dust property profiles 
(from thermal emission mapping), and embedded source distributions (from in- 
frared imaging) we had been seeing in the talks covered the same region of 
space. In other words, since people usually observed their "favorite" object(s) 
with their favorite technique, we almost never knew then density, temperature, 
and velocity field of the same exact set of objects. 1 

During that same Santa Cruz meeting, we calculated that a "Coordinated 
Molecular Probe Line Extinction Thermal Emission" (COMPLETE) Survey 
covering tens of square degrees of sky would take only many weeks of observing 
time, rather than the preposterous many decades of time it would have taken 15 
years earlier. It was obvious that if we could create such a database, and make it 
publicly available, then many researchers (including us!) could finally combine 
all of the powerful techniques developed over the past twenty years, in order to 
measure the physical properties of a carefully chosen set of star-forming regions. 
It was also obvious that we should chose regions that were then about-to-be, and 
now are in the midst of being, observed by the Spitzer Space Telescope under 
the "Cores-to-Disks" (c2d) Legacy Program. 2 By choosing the three Northern 



Elizabeth Lada's pioneering thesis work (ILadal Il99l ILada et all Il99l ILada. Bally, fc Starkl 
I1991TI provides a notable exception to this trend. Through years of work, Lada and colleagues 
mapped the actively star-forming cloud L1630 (a.k.a. Orion B) in both molecular lines (CS, 
with 2' resolution) and with near-infrared cameras (reaching vtik < 13 mag). Lada et al.'s 
work provided the first evidence that massive stars form in clusters, and it also showed that 
the mass spectrum of the gaseous material (self-gravitating or not) is shallower than that of 
stars. COMPLETE, which was not feasible a decade ago, should do for low-mass (fainter) 
star-forming regions what Lada's work did for (brighter) massive star-forming regions, and 
more. For comparison, the total areal coverage of COMPLETE is ~ 20 square degrees, which 
is an order of magnitude larger than the area Lada et al. studied in Orion. 

2 More information about c2d and about COMPLETE (including online data) are available 
through the COMPLETE web site, at cfa-www.harvard.edu/COMPLETE. 
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Hemisphere c2d extended cloud targets (Perseus, Ophiuchus and Serpens), we 
would be assured that a full census of the embedded stellar population and 
its properties would be available around the same time that COMPLETE was 
finished. 

Today, the "COMPLETE" Survey conceived at the infamous Santa Cruz 
coffee break is a healthy two-year old large international collaboration. 1 In the 
winter of 2003/4, we have reached the halfway point in COMPLETE's observa- 
tions, and Spitzer is up and about to begin pouring out the c2d Legacy data. 



2. COMPLETE Observations 

COMPLETE was designed to be executed in two phases. Phase 1, from 2002- 
2004 focuses on observing the "larger context" of the star formation process (on 
the 0.1 to 10 pc scale). All of the Perseus, Ophiuchus and Serpens cloud area 
to be observed with Spitzer under c2d, and slightly more, will be covered by 
COMPLETE in an unbiased way in Phase 1. Phase 2, from 2004-2006, is a 
statistical study of the small-scale picture of star formation, aimed at assessing 
the meaning of the variety of physical conditions observed in star forming cores 
(on the < 0.1 pc scale). In Phase 2, targeted source lists based on the Phase 1 
data are being used, as it is (still) not feasible to cover every dense star-forming 
peak at high-resolution. 

Phase 1 observations (full coverage of all three extended clouds) include: 
FCRAO mapping (40", 0.06 km s -1 resolution) of the molecular gas ( 12 CO and 
13 CO); 2MAS S-based extinction mappi ng (~ 3' resolution) using the NICER 
technique (see lLombardi fe Alvesl . l200ll . and references therein) ; re-processing 
of the 60 and 100 fim ISSA plates to create new maps of dust column density 
and temperature; and SCUBA/JCMT mapping at 850 fim (14" resolution) of 
all regions where Ay > 5. To give a feel for the enormity of this data set, we 
can point out that just one of our molecular line maps of Perseus has ~ 10 5 
spectra of 10 3 channels each. The entire data set is freely available at cfa- 
www.harvard.edu/COMPLETE. 

Phase 2, which has just begun, includes IRAM 30-m spectral line and con- 
tinuum (MAMBO) mapping of selected cores (resolution ~ 10—20"), and higher- 
resolution extinction mapping using both the c2d data itself and ground-based 
near-IR data from 4- to 8-m class telescopes. 



3. Warm Bubbles in "Cold" Clouds? 

Figure 1 shows COMPLETE's velocity-integrated 13 CO map 3 of Perseus overlaid 
on a new dust column density map of the region created using 60 and 100 /im 
ISSA data 4 . The molecular line map is dominat ed by the long chain of clouds 
first mapped out bv lBachiller fc Cernichard ( 19861 "). while the ISSA (IRAS) based 



3 Figure 1 is identical to an image shown at the meeting summarized in this volume. Current 
maps have more coverage, and are available at cfa-www.harvard.edu/COMPLETE. 

4 See lArce fc Goodman! ltl99ll and lRidee et~aT] J2004afl for methodology. 
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map clearly shows a big bubble. 5 Figure 2 shows that the ring dominating the 
60/100-/im-based column density map of Perseus is apparently heated by even 
warmer material interior to it. 

This kind of morphological dissimilarity between warm dust and cold gas is 
also obvious in the famous IRAS image of the Orion molec ular clouds, which i s 
dominated by large shells emanating from evolved stars fe.g. lZhang et all 11989). 
In low-mass star forming regions like Taurus, however, an IRAS-based map is 
nearly indistinguisha ble from an integrated 13 CO map or an extinction map 
((Abergel et allll9 94). In regions exhibiting a mixture of warm and cool dust, 
a 60/100-^m-derived column density map can bear a strong resemblance to the 
cool clouds' morphology, but the dust temperature map often highlights shells. 
For example, our COMPLETE ISSA-based dust temperature map of Ophiuhcus 
(available online) shows a 2-pc diameter heated rin g ~ around p-O ph, centered 
a degree north of the so-called "/9-Oph" cluster (see lLi et alJ . l2004| )). 

In general, how similar an extinction or molecular-line map looks to a dust 
emission map depends sensitively on the mixture of dust temperature and com- 
position along the line of sight. At a given wavelength, the same amount of 
dust at a higher temperature will always emit more than it would at a lower 
temperature. Thus, once significant amounts (> amount of cold dust) of heated 
dust are present in any region, the appearance of that region will be dominated 
by the distribution of that dust in thermal emission maps. Changes in dust com- 
position or size distribution that are correla ted with dust heatin g complicates 
the picture even further and are discussed in lRidge et alJ (|2004ah . 



3.1. Impact Assesment 

How much do warm bubbles impact the star formation process overall? In 
Perseus, it is not even clear that the bub ble and the molecular clouds lie at 
the same distance (see iRidge et al. On the other hand, in Ophiuchus, 

there is clear evidence for interaction of the warm bubble around p-Oph with 
the molecular gas, and it even appears that this bubble may have triggered the 
formation of main star-forming cluster there (confusingly known as the " p-Oph" 
cluster but actually a degree to the South of the star, see lLi et all 12004). 

Simple calculations can show that the energy associated with the shells in 
Perseus and Ophiuchus is enough for either of them to have been generated by a 
~ 10 51 erg supernova going off in dense gas, but our initial investigations favor 
their creation by persistent winds from evolved stars. Whatever their origin, 
though, the energy in these shells may well contribute very substantially to the 
driv ing of turbulence in the den se ISM, as many theoretical simulations demand 
(see iMac Low h Klessenl . 120041 . and references therein) . 

When our Phase 1 molecular-line maps of Perseus (2004) and Ophiuchus 
(2005) are finished, we will measure the kinematic impact on the cold gas of the 
shells we see in warm dust directly, by analyzing the temperature, density and 
velocity structure of the regions most likely to be interacting. Figure 3 shows 
correlations amongst 2MASS(NICER)-, IRAS(ISSA)-, and 13 CO-derived mea- 
sures of column density in Perseus, at positions coincident with the warm ring 



5 The 2MASS/NICER extinction map, available at the COMPLETE web site, looks morpholog- 
ically like the integrated 13 CO map. 
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(white dots) and with the surrounding cold gas (black dots; see Figures 1 and 
2). Notice, in Figure 3b, what a poor tracer of cold material the ISSA map is. 
Notice, in Figure 3a, how W( 13 CO) (oc gas column density) is systematically 
higher in the cold clouds than in the ring for any given 2MASS-derived column 
density, even though the N 2M ASS and PF( 13 CO) are well-correlated overall. In 
order to correctly assess the physical influence of bubble-creating winds on the 
evolution of star formation in molecular clouds, this kind of confounding, con- 
fusing but "Coordinated" Molecular Probe Line Extinction Thermal Emission 
information is, alas, all required. It is not possible to extract relevant measures 
of the mass, temperature, and velocity of the material involved in shell/cloud 
interactions from any one technique alone. 



4. Bipolar Outflows 



The long-term influence of bipolar outflows from young stars on star-forming re- 
gions depends on the spatial and temporal frequ ency of their o ccurrence and on 
the to tal amount (|Margulis. Lada. Snelill988h and variation ( Arce Sz Goodmanl 
2001) of outflow energy with time. We seek to measure all of these parameters 
for the three regions under study in COMPLETE. 

In 2003, we undertook a systematic search for outflows in Perseus, using an 
early version of the huge (> 150,000 spectr a), Nyquist sampled, 1 2 CO map be- 
ing cre ated under Phase 1 of COMPLETE (jFallscheer et all 120031 : iRidge et all 



2004 bl) . By using the Spectral Correlation Function (hereafter SCF. lRosolowskv etall 
1999: 1Padoan~ . joodman, fc Juvelal . l2003h . to identify extended patches of either 
very high or very low spectral correlation in the 12 CO map, we can recover all 
ofl the known outflows in Perseus, as well as many more as "candidate" flows. 
Following up with customized channel maps of the candidate regions has allowed 
us, to date, to roughly double the number of known young stellar outflows in 
Perseus. 

The task that stands before us now is use the power COMPLETE gives us 
to measure density, velocity, and temperature independently (see ^3.1.|) in order 
to correctly determine the energetic impact of bipolar flows. Unlike the case with 
spherical shells, where we would need nearly a galaxy's worth to do a proper 
statistical job, a cloud as big as one of our COMPLETE targets (e.g. Perseus) 
is probably large enough to assess the long-term statistical impact of the bipolar 
flows. "AH" we need to do is to estimate the flows' sculpting and turbulence- 
driving potential now and in the past! We will estimate outflow-driving source 
ages using Spectr al Energy Distributi ons (based on combining ROSAT, 2MASS, 
c2d, BOLOCAM (iKuoch et all I2003T ). and COMPLETE/SCUBA data), allow- 
ing us to study outflow behavior as a function of stellar age. Finally, we will 
construct Monte Carlo models of outflow-cloud interactions where the tempo- 
ral and spatial distribution of outflows, and the density/temperature/velocity 
distribution of the ambient material are all drawn from observations of a single 
star- forming region (e.g. Perseus). 
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5. Dense Cores 

Several groups working on numerical simulations of turbulence in the ISM have 
pointed out that often what looks like a "cloud" or "core" in 2D is really just a 
chance column-density peak, created in 3D by the chanc e superposition of many 
phys ically unrelated structures along the line of sight (c.f. lOstriker. Stone, Gammid . 

In scouring the COMPLETE maps of Perseus in our search for outflows 
(see i!4.|) . the SCF turned up several regions which appear as column density 
peaks, but are in fact comprised of many different velocity components along 
the line of sight. COMPLETE's high velocity resolution allows us to separate 
these velocity features where they would have been blended before. 

Once Phase 1 of COMPLETE is done, we will calculate the fraction of 
column density peaks attributable to these chance coincidences (based on their 
spectra), rather than real 3D density enhancements. In a related effort, we 
will study the spatial dependence of line width on position for the "real" (3D) 
cores and for the "fake" (2D projected) cores, in order to measure how com- 
mon "coherence" (pla teaus of near-thermal line width) is in each group (see 
iGoodman et~aH . 11998ft . We will also compare the distribution of "clumps" be- 
ing found with SCUBA at 850 (described in Johnstone et al. 2004) with the 
cores identified in the spectral-line analysis. 

Using Phase 1 and 2 COMPLETE data, a paper by graduate student Scott 
Schnee et al. will address the context of core "rotation," by analyzing veloc- 
ity gradients and spectral correla t ion in and around dense cores, rather than 
just in the cores (|Goodman et all 11993ft . Paola Caselli and Mario Tafalla will 



lead an inves tigatio n of whether chemical depletion (e.g. ICaselli et all Il999l : 
iTafalla et all l200l 12004ft can be used clS Si clock in dense cores. The high- 
resolution NICER-based infrared extinction maps currently being proposed and 
constructed by Joao Alves et al. will be used as the unbiased tracer of column 
density in many of these core studies. 



6. What if stars move far enough from home to matter? 

Recently, in work not related to COMPLETE, we have determined that PV 
Ceph, a Herbig Ae/Be star, is currently moving at ~ 20 km s^ 1 , and was 
likely ejected, 500,000 years ago, from a clus ter (NGC7023) 10 pc aw ay from 
the molecular cloud where it finds itself now (|Goodman fc Arcel . 12004ft . Other 



researchers ha ve found that stars (still) in clusters commonly w hiz around at 10's 
ofkms -1 fe.g. lPlambeck et all 119951 : feallv fc Beinurth Ll200lft . And, numerical 
simulations predict that anytime stars form in groups-even small groups-some 
are ejected at speeds far in excess of the background velocity d istribution of 
the gas (e.g. iKiseleva et al. |, 119981 : iBate. Bonnell. fc Bromml |2002) . Therefore, 
it is reasonable to worry that stars might move so far beyond the confines of 
the core where they first formed that it would be hard to find all young stars' 
true birthplaces. (For reference, at just 2 km s _1 relative to the gas, a star will 
move from the center to the edge of an 0.1 pc-diameter core in just 250,000 
years.) It is possible that cores move along with the stars they form, relative to 
lower-density material, making star-birthplace pairing less of a problem, but it 
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is doubtful that we will ever be able to associate every young star we see with 
its exact birthplace. 

Even the mild motion of stars with respect to their ancestral homes will 
confound simplistic calculations of star formation efficiency. So we 6 have begun 
a study of star-gas motion aimed at assessing the importance of "stellar scram- 
bling" in a statistical way. Ultimately, we will use COMPLETE data to map 
out the density and velocity distribution of gas, and existing and new radial ve- 
locity and proper motion data to measure stellar velocities. To estimate stellar 
ages, we will use the same kind of SED analysis discussed in When the 
scrambling study is finished, we will be able to assess how far off an estimate 
of "star formation efficiency" based on dividing the mass in young stars by the 
mass of star-forming gas in a particular region on the sky is likely to be. 

7. COMPLETEing the Work of the Three Wise Men 

The Coordinated Molecular Probe Line Extinction Thermal Emission (COM- 
PLETE) Survey of Star-Forming Regions was begun two years before the three 
wise men (Hollenbach, McKee & Shu) turned 60, and it was explicitly designed 
to answer many of the questions those wise men's work have raised but not 
answered. Perhaps a seventieth-birthday party gift will include some answers. 
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